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We explore the consequences of the neutrino mass matrix having a hidden Z2 symmetry and one 
zero eigenvalue. When implemented, these two conditions give relations among the mixing angles. 
In addition, fitting these relations to the existing oscillation data allows limits to be placed on the 
parameter of the symmetry. 
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1. INTRODUCTION 

Neutrino physics can anticipate an era of higher precision measurements with the upcoming generation of neutrino 
experiments. In the past, measurements have shown that the mixing pattern of lepton sector is quite different from 
that of quark sector. In the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) parameterization [l| for lepton mixing there 
are two large mixing angles. The atmospheric mixing angle 9 a = #23 is almost maximal while the solar mixing angle 
S = 612 is also large and the reactor mixing angle 9 X = #31 nearly vanishes. The recent results are summarized in 
Table|H We can see that the uncertainties in mixing angles are not particularly small. In most measurements there 
is roughly a 3° deviation at la confidence level. The mixing matrix which incorporates these angles and diagonalizcs 





A s (l(T 5 eV 2 ) 


A a (HT 3 eV 2 ) 


sin 2 e s (6s) 


sin 2 6a (6a) 


sin 2 9 X (0 X ) 


Central Value 


7.67 


2.39 


0.312 (34.0°) 


0.466 (43.0°) 


0.016 (7.3°) 


1 a Range 


7.48 - 7.83 


2.31 - 2.50 


0.294 - 0.331 
(32.8 - 35.1°) 


0.408 - 0.539 
(39.7 - 47.2°) 


0.006 - 0.026 
(4.4-9.3°) 



TABLE I: The global 3v fit [| for the neutrino mass-squared differences and mixing angles including the available data from 
solar, atmospheric, reactor (KamLAND and Chooz) and accelerator (K2K and MINOS) experiments (j§] . 

the mass matrix M v via U T M V U = M^ lag is given by 



U„ 



-a C 5 S a Sa;6 D C s C a ~h S s S a S x G ' L 




(i) 



with (s a ,c a ) 



(sin# Q , cos#q,) for a = s,a,x. 8d is the Dirac phase and we have neglected Majorana phases. 
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From Table I we see that a good first approximation is to take 9 X = 0, 9 a = 45° which gives 



/ cos( 



U V {0 S ) = 



V 



sin0 s 



^ sin( 



sin u s 
\ cos 9 S 

\ cos 9 S 



\ 



Vl V 2 V 3 



(2) 



Using this as a starting point, we wish to investigate whether there is an underlying symmetry G of the neutrino 
mass matrix. This matrix must satisfy [G, M„] = 0, or G T M V G = M v . Given a G, the transformation GU U also 
diagonalizes M v . But U v is unique except for phases. This can be seen by supposing that d v is a unitary matrix 
such that U v d v also diagonalizes M v . For this to be true, d v must satisfy d^ M^ ms d„ = M^ lag and this implies that 
= 1. Since GU U diagonalizes M v , it must have the form GU V = U v d vi where d v has 1 or —1 diagonal elements. 
Thus 



UlGU v = d v 



\ 



& G = U v d„Ul . 



(3) 



There are only eight possible combinations for the elements of d v . Two of these are the unit matrix and its negative, 
both of which define G as a multiple of the identity. Of the remaining six, three have two entries of +1 and one of 
— f , while the other three have two entries of —1 and one of +1. These diagonal matrices differ by an overall minus 
sign, so only one of the two types is independent. If we choose the three with one entry of 1 and two entries of —1 
(det(G) = 1), then it is easy to see that multiplying any pair of these diagonal matrices will result in the remaining 
matrix. Hence, in reality, only two of these matrices are independent and both represent a Z2 symmetry Since the 
independent G's commute, the horizontal symmetry of lepton mixing is Z2 x Zi if neutrinos are Majorana fermions 

m- 

A representation of G can be obtained using 



G = dxviv{ + d 2 v 2 vl + d 3 v 3 vl . 
Since the eigenvalue 1 can occur in three places, there are three symmetry matrices G 



G 2 = 



(4) 



and 





V2s s c s 


y/2s s c s 


G\ = V2s s c s 




si 


\ \/2s s c s 


S 2 


-cl 






G 3 = 





where the subscript i on Gi denotes the component of dv that is +1. G3 gives /x — t symmetry 0] while G± is 
symmetric and commutes with G3. For simplicity we can parameterize the solar mixing angle 9 S as 

-k . „ y/2 



sin 9* 



Then the mixing matrix ^ takes the form [?J 



U v {k) 



V2+k' 2 
1 



-y/2 
-k 



Vfc 2 



\ 



(5) 



VWk? ^2(2+fe 2 ) v/2 
1 -fc 1 

^2(2+fe 2 ) V2 J 



U k . 



(6) 
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Consequently, the symmetry transformation matrix Gi(9 s ) can be reexpressed in terms of k 



2 k 2,h 2.L 

k 2 



Gi(k) = | 



Although we can "derive" a generalized form of G\ symmetry transformation matrix ([7]) given the mixing matrix 
([5]), this relationship cannot be reversed. The mixing matrix U v can not be uniquely determined solely by G\ due to 
the fact that G\ has degenerate eigenvalues. 

Invariancc under Gz requires 6 X = 0°, 6 a = 45°, but invariance under G\ does not, so below we assume that 
the neutrino mass matrix is invariant under G\, not only in the approximation X = 0°, 9 a = 45°, but for general 
values of all the mixing angles. In the next section we use this assumption in the form of Eq.([7]) with general values 
of k to derive relations among the mixing angles. In Sec. 3 we compare our results with the experimental values and 
in Sec. 4 we summarize. 



2. INVARIANCE UNDER THE Z 2 SYMMETRY Gi 



In this section we show explicitly the consequences of generalized G\ symmetry. Only two mass square differences 
have been measured and the neutrino's mass scale has not been determined by experiments. It is possible that one of 
the mass eigenvalues vanishes. This is also theoretically motivated by minimal seesaw model 0. We will explore the 
joint consequences of one vanishing mass eigenvalue and G\ invariance. For simplicity we will postpone discussion 
of CP phases to a later article ■ 



2.1. Constraints on Mass Matrix Elements 

If the neutrinos are Majorana fermions, their mass matrix must be symmetric. We will consider the case that 
there are three generations of light neutrinos. Then, the most general form of the neutrino mass matrix can be 
parameterized as 





(a 


B k 


B 2 \ 


M v = 


B x 


C x 






\b 2 


D 


c 2 J 



(8) 



which has six independent matrix elements. We assume M v is invariant under the G\ symmetry transformation, 



With the help of © and 



G\M V G X = M v . 

Eq. ((HJ) gives two conditions on the neutrino mass matrix elements of 

gi +B 2 k 

Ci+C 2 + 2D-2A ~ k 2 - 2 ' 
B 1 -B 2 1 
C\ - C 2 k 
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(9) 

(10a) 
(10b) 



2.2. Eigenvalues and eigenstates 



If there is a vanishing mass eigenvalue m, = the corresponding mass eigenstate, which can be denoted as 
v = (a,/3,7) T , must satisfy 



4 



(A 


Bi 


B 2 




Ci 


D 


\B 2 


D 


c 2 




(11) 



If we assume a ^ then we get three equations 



A = -pB 1 -aB 2 , (12a) 
Bi = -pCi-aD, (12b) 
B 2 = - P D-aC 2 , (12c) 

where p = /3/a, a = j/a. Thus we have two sets of conditions, (jXTJJ) from G± invariance, and (|12[) from the vanishing 
mass eigenvalue. 

From the relations (fl"2"|) we can express the matrix element A in terms of C\ , C 2 and D 

A = p 2 C 1 +a 2 C 2 + 2paD. (13) 

Now let us use these in the Z 2 relations. Eq. QlOap and Eq. (|10b[) give 

(cr fc + 1)C 2 - (p fc + l)Ci + - o-)Z> = , (14a) 
{pk + l)(2p - k)d + (cr k + I) (2a - k)C 2 + [(2 - k 2 )(p + a) - 2k(l - 2pa)]D = . (14b) 

The above two relations can be reexpressed in terms of only two matrix elements, D and C 2 or C\ respectively 

(p + <j-k)[(crk + l)C 2 + (pk + l)D} = 0, (15a) 
(p + a -k)[(pk + l)C 1 + (ak + l)D] = 0. (15b) 

The mass eigenvalues that are nonzero are given by 

1 



m ± = - 



A + d + C 2 ± y/(A + d + C 2 f + 4(^2 + a 2 + 1) (D 2 - dC 2 ) , (16) 



where we have used (|12bj) and ()12cp and (fT3|) . From (|15al) and (j!5bj) it is obvious that one possible solution to the 
equations for C\,C 2 ,D is 

d = ~^jD, (17a) 
pk + 1 

C 2 = -^±1^. (17b) 
cr fc + 1 

This makes D 2 = C\C 2 and consequently m_ given above would also be zero. Since the experimental data shows 
that two mass square differences between the three neutrino mass eigenvalues are nonzero, we need at least two 
masses nonzero in order to have two oscillation lengths. 

A second solution of (|15a|) . (|15bl) is p = a = — 1/fc but then the three relations in (fT2"j) simply reproduce the 
conditions (JTUJ). So the conclusion is that we must have 

p = k-a (18) 

and the conditions (|14a|) or (|14b[) reduce to an equation for a 

(1 + fc 2 )Ci - C 2 - k 2 D 



k(d +d- 2D) 



(19) 



This relation represents the constraint from G\ invariance which was originally expressed as (|10p where there were two 
independent relations. Using (|18|) . these two relations are satisfied simultaneously and reduce to a single constraint 

m 



The condition (|18p can also be substituted into (|12a[) , (|12bp , and (|12cp to give 

A + kBx . , 

(20a) 

(20b) 

(20c) 

respectively. These three relations are a manifestation of vanishing mass eigenvalue. We can set these equations for 
a equal to get relations among the matrix elements A, . . . , D in terms of the parameter k. Not all of these equations 
are independent but two different relations are possible: 

(A + kB 1 )(C 1 -D)-{B 1 +kC 1 )(B 1 -B 2 ) = 0, (21a) 
(A + kB 1 )(D-C 2 )-{B 2 +kD)(B 1 -B 2 ) = 0. (21b) 





B 2 ' 


Bi + 


fed 


Ci - 


- D 


B 2 + 


kD 


D — 


c 2 



In the next subsection we will write A, . . . , D in terms of the mixing angles and thereby get two relations among 
the mixing angles, again involving k. 



2.3. Reconstruction of Neutrino Mass Matrix 



Using U v from Eq. in M v = U*M^W and comparing with © we get P3 



A = 


2 2 ,2 2 ,2 

c x c s mi+ c x s s m 2 + s x m 3 




(22a) 


Si = 


C-x\$ sCsC-a ^x^aC s \fTli C x [s s C s C a H~ S x S a S s 


m 2 + c x s x s a m 3 


(22b) 


B 2 = 


C-x\s s^s^a H~ ^a;C a C s ]?7li C x [s s C s S a S x C a S s 


\m 2 - s x c x c a m 3 


(22c) 


Cx = 


(s s c a - s x c s s a ) 2 mi + (c s c a + s x s s s a ) 2 m 2 4 


■ c ls 2 a m 3 


(22d) 


c 2 = 


(s s s a + s x c s c a ) 2 m 1 + (c s s a - s x s s c a ) 2 m 2 4 


c 2 x c 2 a m 3 


(22e) 


D = 


($sSa S x CsCa) (^s^a Sx^s $a )^1 ~t~ (Cs-5 a 


SxSsC a )(c s C a + s x s s s a )m 2 - c x s a c a m 3 


(22f) 



where, as mentioned above, we have deferred consideration of CP violation to a later article. The mass eigenvalues 
can be further parameterized in terms of experimentally measured mass square differences: mi = too, m 2 = moVT+r 
and to 3 = for inverted mass hierarchy and mi = 0, m 2 = rriQ^fr and m 3 = mo for normal mass hierarchy where 
mo = V A a and r = A s /A a which is positive. 



2.4. Correlations between Mixing Angles 



To get relations between mixing angles we can substitute (|22p into (|2"Tj) which gives 

- c a c x [c x (c a - s a ) + ks x ] mim 2 = 0, (23a) 
s a c x [c x (c - 8 a ) + ks x ]mim 2 = 0. (23b) 

where we have assumed the mass hierarchy is inverted with vanishing m 3 and nonzero mi and TO2, while c a , s a and 
c x are also nonzero. The only possible solution is 



G 



where the last factor comes from expanding the two mixing angles 9 a and 9 X around the approximations 45° and 0°, 

a = J + S a , 9 X = S X . (25) 

With this solution for k and the reconstructed mass matrix elements (|2"2l substituted back into (|20a|) we find 

« -V2S X . (26) 

Since 6 X is quite small, a should be very large according to (|26[) . 

We still have the condition from Z 2 , Eq. ([19]). Together with ([26]) and <[22j) as well as ([24]) it gives 



(j k c a c x c a 



4 - (2 + 2sl)c a s a f c a + s a Y i f c, - 5 n 1 N 

I c a - s a I c„ + s n s s 



There are two possible solutions 



tan (9 S 



tan#, = s T = — 



^(^a ~t~ ^a) ^a; 
Qz:(Ca "T" ^a) 



(28) 



(29) 



These relations between mixing angles can be used to predict the not well measured X in terms of the solar and 
atmospheric mixing angles. For example (|29[) gives 



C-a Sa Ss 
£ a T Ca C-s 



tan0,d n 



(30) 



Since 8 X is the focus of next generation of neutrino experiments, we use ([30]) to estimate its value. The scatter plot 
is shown in Fig. [TJ A scatter plot based on (|28p would look similar with a steeper slope for the points. 




46 48 



FIG. 1: Prediction of 8 X in terms of 9 a and 9 3 at the 90% C.L. The vertical solid line denotes the experimentally measured 
central value of the atmospheric mixing angle 9 a . 



Another way of expressing the results is to write all of the mixing angles in terms of the parameters a and k. 
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Using z = l/cr, this gives 

z 2 

sin 2 9 X = —z ——: — (31b) 

(k 2 + l)z 2 - 2kz + 2 v ; 

£; 4 z 2 -2/fc 3 z + 2/fc 2 (z 2 + l)-4fcz + 4 1 ' 

sin 2 - (fc 2 + l)z 2 -2fcz + 2 

s ~ fc 4 z 2 -2fc3z + 2fc 2 (z 2 + l)-4fcz + 4' 1 j 

Note that these equations are all unchanged under k, z — > —k, —z so only the absolute value of k can be determined. 



and either 



or 



3. FIT TO EXISTING DATA 



The solutions (|28[l or (f29|) . which give (|31cj) or (|31d|) . arc identical in the following sense - oscillation experiments 
measure sin 2 29 and thus can't distinguish between 9 and 7r/2 — 9. Further, tan(7r/2 — 9) = 1/ tan#, so a fit with 
(|28|) and 9 S assumed greater than ir/4 is identical to one with (|29|) and 9 S assumed less than 7r/4. Having noted this 
we will proceed to fit both (|31cp and (|31dj) with 9 S < ir/4. 

Using Eqs. (|3Ta| . ([Mb} , and (f31c|) . the fit to the data from Ref. gives x 2 mia = 2 - 10 ; l^lmin = 2.09 and z = 0.066. 
At the minimum values of \k\ and z, sin 2 (6» a ) = 0.426 (6 a = 40.7°), sin 2 (9 x ) = 0.0025 (9 X = 2.87°) and sin 2 (9 s ) = 
0.313 (6» s = 34.0°). The 68.3% and 90% confidence contours are shown in Fig.[2] 



0.2 




1.9 



2.1 
k 



2.3 



FIG. 2: The 68.3% and 90.0% confidence contours for the fit using Eqs. (|31a[l . (|31b[l . and (|31c[) are shown in red and orange, 
respectively. The (black) dot indicates the x 2 minimum. 



The distributions of this set of mixing angles arc obtained from the likelihood distribution 

i e -(x 2 (M)-xL)/2 i (32) 

where A is a normalization constant, using 

dP 



dsin 2 (9) 



dk dz 5{sm 2 (9) - f(k, z))Ae~ {x (M)-*™)/ 2 ; (33) 
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where f(k, z) is one of the functions on the righthand side of Eqs. ([31]) . The results are shown in Figs. ([3]). 




FIG. 3: The distributions of the sin 2 Oi obtained using Eqs. (|31a[) , (|31b|) , and (|31c|l are shown. 

As would be expected, none of the distributions is exactly Gaussian. The largest contribution to the minimum % 2 
is associated with sin 2 (# a ) and the influence of terms beyond the quadratic expansion of x 2 (fc, z) can be seen in the 
shape of this distribution. 

If we use Eqs. (|31a| . (|31b| and (|31d|) . the fit to the data has two local minima. The lowest of these gives Xmin = 
0.506, |fc| min = 0.942 and z min = 0.152. At the minimum values of |fc| and z, sm 2 (9 a ) = 0.423(40.5°), sin 2 ^) = 
0.013(6.55°) and sin 2 (6» s ) = 0.311 (33.9°). At the other minimum, where x 2 = 2.73, sin 2 6 S and sin 2 6 X are slightly 
different, but sin 2 9 a = 0.567. This is reflected in the individual mixing angle distributions. The confidence contours 
for this case are shown in Fig. 01 
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FIG. 4: The 68.3% and 90.0% confidence contours for the fit using Eqs. (|31a|) . (|31b[) and (|31d|l are shown in red and orange, 
respectively. The (black) dot indicates the \ 2 minimum. 

The distributions of this set of mixing angles are shown in Figs. ([5]). Here, too, the largest contribution to the 
minimum % 2 is associated with sin 2 (6* a ) and the effect of the second local minimum this is reflected in the distortion 
on the high side of the probability distribution. 

Alternately we can fit for k using the values of sin 2 9 a and sin 2 9 S from Ref . [H but replace sin 2 9 X with the value for 
sin 2 (26» a ) sm 2 (29 x ) published by the MINOS collaboration [l]|. They report sin 2 (26» a ) sin 2 (20 x ) ~ 0.18 ± 0.13 for 
inverted hierarchy and, for normal hierarchy, ~ 0.11 ±0.09. From the inverted hierarchy result we get \k\ = 2.10±0.10 
with a x 2 of 1.86 or |fc| = 0.94 ± 0.15 with a X 2 of 1.35. 

This was all for inverted hierarchy Normal hierarchy, with ml equal to zero, gives, after a lot of work, exactly the 
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0.35 0.40 0.45 0.50 0.55 0.60 0.65" 



sin" 8 S 



dsn 


2 9* 


40 




M\ 




/" 




/ 10 





-0.02 - 0.01 0.00 0.01 0.02 0.03 0.04 0.05 



FIG. 5: The distributions of the sin 9i obtained using Eqs. (|31a[) , (|31b|) and ()31d|) are shown 



results of inverted hierarchy, (|24[) . (|28|) . (|29|) . The parameter er is a different function than (|26|) . 



s a (l + k 2 ) - c a 

fc(So + C a ) 



(34) 



but this just amounts to a reparameterization of ([31]) with no physical consequence. Using the MINOS number for 
normal hierarchy we find the same values, including the errors, for as for the inverted hierarchy MINOS number. 
The x 2 values are smaller at 1.42 or 0.80. 

With either MINOS value and for either value of |fc| the fitted value of sin 2 9 S is stable at 0.312, the fitted value 
of sin 2 9 a varies only slightly from 0.46 for the larger |fc| to 0.42 for the smaller value, but sin 2 8 X is less than 0.001 
for the larger \k\ but equal to 0.015 for the smaller. 



4. SUMMARY 



A hidden Zi symmetry, as given by Eq. ([5]), results in only two possible sets of conditions on the neutrino mixing 
angles. Assuming 8 S < -k/A then either 



tanfl, 



with confidence contours shown in Fig. [2] or 



, {Sa c a) i 



(35a) 
(35b) 



tan 0, = 



k 



(36a) 
(36b) 



with the confidence contours shown in Fig. [4] 
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